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Abstract The strength and effect of Coulomb coiTelations in the (superconducting when 
hydrated) x ~l/3 regime of Naa;Co02 have been evaluated using the correlated 
band theory LDA+U method. Our results, neglecting quantum fluctuations, are: 

(1) there is a critical Uc = 3 eV, above which charge ordering occurs at x=l/3, 

(2) in this charge-ordered state, antiferromagnetic coupling is favored over ferro- 
magnetic, while below Uc, ferromagnetism is favored; and (3) carrier conduction 
behavior should be very asymmetric for dopings away from a::=l/3. For x < ^, 
correlated hopping of parallel spin pairs is favored, suggesting a triplet super- 
conducting phase. 



1. Introduction 

Since the discovery of high temperature superconductivity in cuprates, there 
has been intense interest in transition metal oxides with strongly layered, (quasi) 
two-dimensional (2D) crystal structures and electronic properties. For several 
years now alkali-metal intercalated layered cobaltates, particularly Na2:Co02 
(NxCO) with X ~ 0.50 — 0.75, have been pursued for their thermoelectric 
properties. [1] Lia;Co02 is of course of great interest and importance due to 
its battery applications. The recent discovery [2] and confirmation[3-5] of su- 
perconductivity in this system, for x 0.3 when intercalated with H2O, has 
heightened interest in the NxCO system. 

The crystal structure [6-8] is based on a 2D C0O2 layer in which edge- 
sharing CoOg octahedra lead to a triangular lattice of Co ions. Na donates its 
electron to the C0O2 layer, hence x controls the doping level of the layer: x=0 
corresponds to Co^+, S=^ low spin ions with one minority t2g hole, and x = 1 
corresponds to non-magnetic Co^+. Nearly all reports of non-stoichiometric 
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materials quote values of x in the 0.3 - 0.75 range, and the materials seem gen- 
erally to show metallic conductivity. Reports of the magnetic behavior are of 
particular interest to us. For x in the 0.5 - 0.75 range, the susceptibility x{T) 
is Curie- Weiss-like (C-W) with reported moment of the order of magnitude 1 
jjLB per Co^+[2, 3] which indicates the presence of correlated electron behavior 
on the Co ions. Magnetic ordering at 22 K with very small ordered moment 
has been reported for x=0.75[9] and Wang et al. measured field dependence 
[4] that indicated the spin entropy of the magnetic Co system is responsible for 
the unusual thermoelectric behavior. Thus for x > 0.5 magnetic Co ions and 
magnetic ordering give evidence of correlated electron behavior. 

However, for H2O intercalated samples with x «0.3, {i.e. the superconduct- 
ing phase) C-W behavior of x vanishes. [3, 5, 10, 1 1]. It is extremely curious 
that the appearance of superconductivity correlates with the disappearance of 
Co moments in the samples. From a single-band strongly interacting view- 
point, the X = system corresponds to the half-filled triangular lattice that has 
been studied extensively for local singlet (resonating valence bond) behavior. 
[12] The X 0.3 region of superconductivity in NxCO is however well away 
from the half-filled system, and the behavior in such systems is expected to 
vary strongly with doping level. 

There is now a serious need to understand the electronic structure of the nor- 
mal state of the unhydrated material, and its dependence on the doping level 
X. The electronic structure of the a;=l/2 ordered compound in the local den- 
sity approximation (LDA) has been described by Singh. [13] Within LDA all 
Co ions are identical ("Co^'^+"), the Co t2g states are crystal-field split (by 
2.5 eV) from the eg states, and the t2g bands are partially filled, consequently 
the system is metallic consistent with the observed conductivity. The t2g band 
complex is 1^ ^ 1.6 eV wide, and is separated from the 5 eV wide O 2p 
bands that lie just below the Co d bands. Singh noted that the expected on-site 
Coulomb repulsion U =5-8 eV on Co gives U » W and correlation effects 
can be anticipated. A value of ?7 Ri 4 eV was assumed by Wang, Lee, and Lee 
[14] to justify a strongly correlated t — J model treatment of this system. While 
it must be kept in mind that the study of this system is still in its infancy and 
no clear experimental data plus theoretical interpretation agreement has estab- 
lished the degree of correlation, we also take the viewpoint here that effects of 
on-site repulsion need to be assessed. 

Although the experimental evidence indicates nonmagnetic Co ion in the 
superconducting material, most theoretical approaches [14-1 6] consider the 
strongly interacting limit where not only is U important, it is large enough 
to prohibit double occupancy, justifying the single band t — J model. Another 
question to address is whether the single band scenario is realistic: indeed the 
rhombohedral symmetry of the Co site splits the t2g states into ag and e'g rep- 
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resentations, but the near-octahedral local symmetry leaves their band centers 
and widths very similar. 

In this paper we begin to address the correlation question using the corre- 
lated band theory LDA-i-U method. We focus on the x « 1/3 regime where su- 
perconductivity emerges. We find that U >Uc=3eV leads to charge ordering 
at x=l/3 accompanied by antiferromagnetic (AFM) spin order; of course, the 
fluctuations neglected in the LDA-i-U method, the availability of three distinct 
sublattices for ordering, and the tendency of the Na ions to order[8] (which can 
mask other forms of ordering at the same wavevector), can account for the lack 
of ordering (or of its observation). 

2. Method of Calculation 

Two all-electron full-potential electronic methods have been used. The 
full-potential linearized augmented-plane-waves (FLAPW) as implemented in 
Wien2k code [17] and its LDA-i-U [18] extension were used. Valence and con- 
duction s,p, and d states were treated using the APW-i-lo scheme [19], while 
the standard LAPW expansion was used for higher Vs. Local orbitals were 
added to describe Co 3d and O 2s and 2p states. The basis size was determined 
by RmtKmax = 7.0. In addition, the full-potential nonorthogonal local-orbital 
minimum-basis scheme (FFLO) of Koepemik and Eschrig[20] was also used 
extensively. Valence orbitals included Na 2s2p2>s2>p'id, Co 3s3p4s4j»3(i, and O 
2s2p3s3p3d. The spatial extension of the basis orbitals, controlled by a confin- 
ing potential [20] (r/ro)*^, was optimized for the paramagnetic band structure 
and held fixed for the magnetic calculations. The Brillouin zone was sampled 
with regular mesh containing 50 irreducible k-points. Both popular forms[21] 
of the LDA-I-U functional have been used to assess possible sensitivity to the 
choice of functional, but in these studies the differences were small. We do not 
consider interlayer coupling in the work presented here, which allows us to use 
a single layer cell in the calculations. 

3. Results of Self-Consistent Calculations 

LDA electronic structure at x = \. The crystal field sphtting of 2.5 eV 
puts the (unoccupied) Cg states (1 eV wide) well out of consideration. The 
trigonal symmetry of the Co site splits the t2g states into one of ag symmetry 
and a doubly degenerate e'g pair. The ag band is 1.5 eV wide (corresponding to 
f = 0.17 eV in a single band picture) while the e'g states have nearly the same 
band center but are only 1.3 eV wide; hence they lie within the ag band. As 
might be anticipated from the local octahedral environment, there is mixing of 
the Qg and e'g bands throughout most of the zone, and the ag DOS does not 
resemble that of an isolated band in a hexagonal lattice. For the paramagnetic 
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Figure 1. Effect of the intraatomic repulsion U on the magnetic moment of the Col and Co2 
ions for ferromagnetic order. The LDA+U method in the FPLO code was used. Disproportion- 
ation to formal charge states Co'^^ and Co*^ states occurs above Uc = 3 eV. 



case X = ^ corresponds to | filling of the t2g band complex, resulting in hole 
doping into the e'g states as well as in the ag states. 

Singh found that ferromagnetic (FM) phases seemed to be energetically fa- 
vored for all noninteger x [13]. No ferromagnetism is seen in this system, so 
NxCO becomes another member in small but growing list of compounds[22] 
whose tendency toward FM is overestimated by LDA because they are near a 
magnetic quantum critical point. We confirm these FM tendencies within LDA 
for x=l/3, obtaining a half metallic FM state with a moment of l/i^/Co that 
is distributed almost evenly on the three Co ions, which occupy two crystallo- 
graphically distinct sites because of the Na position. The exchange splitting of 
the t2g states is 1.5 eV, and the Fermi level (Ep) lies within the metallic mi- 
nority bands and just above the top of the fully occupied majority bands. The 
FM energy gain is about 45 meV/Co. With the majority bands filled, the filling 
of the minority t2g bands becomes |, leading to larger e'g hole occupation than 
for the paramagnetic phase. We conclude that, in opposition to much of the 
theoretical speculation so far, x = | is necessarily a multiband {ag + e'g) sys- 
tem. Our attempts using LDA to obtain self-consistent charge-ordered states, 
or AFM spin ordering, always converged to the FM or nonmagnetic solution. 

LDA+U Magnetic Structure and energies First the behavior of the LDA-i-U 
results versus U were studied (on-site exchange was kept fixed at 1 eV). The 
dependence of the magnetic moment on U (obtained from the FPLO code) for 
FM ordering is shown in Fig. 1. Recall that the ordered array of Na ions in 
our cell gives two crystallographically inequivalent Co sites. For U < Uc ^ 



Figure 2. <Upper left> The charge ordered triangular Co lattice honeycomb lattice, with 
antiferromagnetic spin order designated by solid circles (t), dashed circles (|), and filled circles 
(Co^+S* = sites). Lines denote nn magnetic couplings. <Upper right> addition of a t electron 
converts a Co^^ S — ^ site to a nonmagnetic site. Hopping of a neighboring hole to this site 
costs 4J in energy. The lower two panels illustrate the correlated pair hopping process after a | 
hole is added to the system: hopping of the hole to a neighboring site, followed by refilling of 
the site by the added hole, results in an identical (but translated) state. 



3 eV, the moments are nearly equal and similar to LDA values (which is the 
U -^0 limit). Above Uc, charge ordering accompanied by a metal-insulator 
(Mott) transition occurs by disproportionation into nonmagnetic Co'^^ and two 
S = ^ Co^+ ions. Nonmagnetic Co'^+ states lie at the bottom of the 1 eV 
wide gap, with the occupied Co^+ e'g states 1-2 eV lower. After adding the 
on-site correlation to LSDA results, the hole on the Co^+ ion occupies the ag 
orbital. A possibility that we have not attempted would be to obtain a solu- 
tion in which the hole occupies an e'g orbital, in which case one should then 
investigate orbital ordering in addition to charge- and spin-ordering. 

Reasonable estimates for cobaltates put U at 5 eV or above, so we now 
concentrate on results for C/=5.4 eV, which we expect is near the lower end 
of reasonable values. For this value of U both FM and AFM spin-ordered 
solutions are readily obtained, with AFM energy 1.2 mRy/Co lower than for 
FM order. In terms of nn coupling on the charge-ordered honeycomb lattice, 
the FM - AFM energy difference corresponds to J = 11 meV. Referring to the 
paramagnetic bandwidth identified above, the corresponding superexchange 
constant would be At^ /U ~ 20 meV. 

Discussion and Comparison with Experiment These calculations establish 
that at x=l/3, there is a strong tendency to charge-order, and that there is a nn J 
of antiferromagnetic sign; hence we consider as reference the \/3 x ^/?, charge- 
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ordered AFM state shown in Fig. 2. Considering the charge-ordering energies 
as dominant over the magnetic energies, the fundamental problem at .x = ^ be- 
comes the spin behavior of the half-filled honeycomb lattice. Spin correlations 
and quantum fluctuations on the honeycomb lattice have been considered by 
Moessner et al. [23] based on the quantum dimer model, where singlet-pairing 
regimes indeed are found in which the rms magnetization on a site is strongly 
reduced. Such pairing would strongly suppress the Curie-Weiss susceptibility. 

The foregoing discussion neglects (among other aspects) the metallic nature 
of NxCO. Now we consider "doping" away from x=l/3. Addition of an elec- 
tron (of, say, spin up) converts a Co^"*" | to a Co^"*", that is, it destroys a spin 
down hole which also was a potential carrier (if charge order were lost). This 
frees up a site for hopping of a neighboring hole, but the energy cost of doing 
so is 4 J (loss of two favorable J and gain of two unfavorable J) and thus is 
strongly inhibited. Now we turn to the removal of an f electron (addition of 
a hole) corresponding to superconducting region x < |, which has quite a 
different effect. This type of doping converts an inert Co^+ to a Co^+ that is 
surrounded by six Co^+ sites with alternating spins. Single particle hopping 
is disallowed (strictly speaking, it costs U); however, correlated parallel-spin 
pair hopping as shown in Fig. 2 returns the system to an equivalent state and 
therefore requires no net energy. This process suggests a tendency toward 
triplet pairing superconductivity in this regime, although our results are too 
preliminary to permit serious conclusions. Although some progress on such 
processes might be made analytically within a single band model, the multi- 
band nature of NxCO turns even the relatively simple Hubbard model on a 
triangular or honeycomb lattice into a formidable numerical problem. 

4. Summary and Acknowledgments 

Now we summarize. We have used the LDA-i-U method to evaluate the ef- 
fects of Hubbard-like interactions in NxCO, and find charge disproportionation 
and a Mott insulating state at a; = ^ when fluctuations are neglected. Nearest 
neighbor coupling J Kill meV provides AFM correlations. Indications based 
on this "mean field" AFM charge-ordered state are for very different behavior 
for electron or for hole doping relative to x = ^ ; hole doping from this point 
tends to favor parallel-spin pair-hopping and thus possible triplet superconduc- 
tivity. Fluctuation effects may however be substantial. 

We acknowledge important communications with R. T. Scalettar, R. R. R 
Singh, R. Cava, B. C. Sales, and D. Mandrus. J. K. was supported by Na- 
tional Science Foundation Grant DMR-0114818. K.-W. L. and W. E. R were 
supported by DOE Grant DE-FG03-01ER45876. 
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